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ABSTRACT 

Gravitational lenses offer the possibility of accurately determining the Hubble pa- 
rameter (Ho) over cosmological distances, and B0218+357 is one of the most promising 
systems for an application of this technique. In particular this system has an accu- 
rately measured time delay (10.5 ± 0.4 d; Biggs et al. 1999) and preliminary mass 
modelling has given a value for H of 69±^ kms" 1 Mpc" 1 . The error on this estimate 
is now dominated by the uncertainty in the mass modelling. As this system contains 
an Einstein ring it should be possible to constrain the model better by imaging the 
ring at high resolution. To achieve this we have combined data from MERLIN and the 
VLA at a frequency of 5 GHz. In particular MERLIN has been used in multi- frequency 
mode in order to improve substantially the aperture coverage of the combined data 
set. The resulting map is the best that has been made of the ring and contains many 
new and interesting features. Efforts are currently und erway to exploit the new da ta 
for lensing constraints using the LensClean algorithm ( Kochanek fc Narayan 1992] ). 



1 INTRODUCTION 



The gravitational lens system B0218+357 (Patnaik et al 



|199S ) has long been recognised as an excellent system for 
measuring the Hubble constant (Ho) over cosmological dis- 
tances using the method of Refsdal (1964). The system has a 
simple morphology comprising two compact images (A and 
B) of a flat-spectrum radio core and a radio Einstein ring 
(Figs [l], ^ and ^). There is also a kpc-scale radio jet that 
extends south from the bottom of the Einstein ring (seen 
most clearly in Fig ^). Both lens and source redshifts are 
well determined , at 0.6847 (Browne et al. 1993) and 0.96 
(Lawrence 1996) respectively. In addition, the time delay 
between the two compact components has bee n accurately 
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A free parameter in the model is the position of the lens- 
ing galaxy centre which is poorly determined from Hubble 
Space Telescope (HST) optical and infra-red observations. 
Lehar et al. (2000) have pointed out that the uncertainty 
in the position of the lensing galaxy implies an uncertainty 
in Hq that is much greater than that given in Biggs et al. 
(1999). 



meas ured (10.5 ± 0.4 d) from VLA monitoring (Biggs et al 
199£). Preliminary modelling of the mass distribution in the 
lensing galaxy usi ng a Singular Isothermal Ellipsoidal (SIE ) 
parameterization ( Kormann, Schneider & Bartelmann 1994 ) 
gave a value for Ho of 69^g kms _i Mpc" 1 . The quoted er- 
ror is a 95 per cent confidence limit (statistical) . Constraints 
on this model come from the VLBI substructure of the com- 
pact images (Patnaik, Porcas & Browne 1995) as well as 
the flux density ratio measured from the VLA monitoring. 
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A potential source of model constraints in B0218+357 is 
the radio brightness distribution in the Einstein ring which 
has been neglected until now. Einstein rings have proved to 
be a valuable source of constraints in other lens systems (e.g. 
MG 1131+0456 and PKS 1830-211) as they probe the mass 
distribution in the lensing galaxy at many points, thus pro- 
viding many more constraints than are available from lenses 
t hat do not cont ain images of large-scale extended structure 
(Kochanek 1990). Techniques that have been developed to 



optimise lens models using radio maps of Einstein ring s in- 
clude the R ing Cycle algorithm ( Kochanek et al. 1989) and 
LensClean (Kochanek & Narayan 1992). 



The Einstein ring in B0218+357 is believed to be an 
image of part of the extended emission of the kpc-scale ra- 
dio jet. It is also the smallest known Einstein ring with a 
diameter of only 335 milliarcsec (mas), equal to the separa- 
tion between the two compact components. Because of its 
small size, producing reliable maps of the brightness distri- 
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bution in the ring that have high resolution and high sen- 
sitivity is difficult. The previous best map was made at a 



freq uenc y of 5 GHz with the MERLIN array ( Patnaik et al i 
199 3j), \ >ut due to the incomplete aperture coverage of the 
observations no attempt was made to interpret details of 
the Einstein ring in a lensing context. This paper presents 
a new map of B0218+357 made from data obtained with 
both MERLIN (in multi-frequency mode) and the VLA at 
frequencies around 5 GHz. The excellent aperture coverage 
when these data are combined has allowed us to produce the 
best map of this system at this frequency to date. Work is 
currently underway to exploit the new map for model con- 
straints (Wucknitz et al., in preparation). 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 MERLIN data 

Observations were taken in June 1995 with the MERLIN ar- 
ray in multi-frequency synth esis (MFS) mode (see Conw ay, 
Cornwell & Wilkinson (]l990[) or Sault & Conway (]l999[) for 



a full discussion of this technique). The advantage of this 
technique is that for a given baseline separation, changing 
the observing frequency changes the spatial frequency of the 
measured visibility and hence the position of the visibility 
in the (u, v) plane. This allows the (u, v) plane to be filled in 
more fully than would be possible for observations taken at a 
single frequency with the same array and covering the same 
period of time. This in turn allows the sky brightness distri- 
bution to be reconstructed more accurately. A complication 
arises from the spectral index of the emission which, if uncor- 
rected for, will produce imaging artefacts and significantly 
degrade the dynamic range (the ratio of peak brightness to 
off-source error) of a map. 

For these observations three frequencies around 5 GHz 
separated by 320 MHz were used. These were 4546, 4866 and 
5186 MHz, a frequency range of ±7 per cent. The bandwidth 
of the individual frequencies was 16 MHz. B0218+357 was 
observed alte rnately with the nea rby JVAS calibrator source 
B0233+359 (Patnaik et al. 1992) for phase-referencing pur- 
poses, a complete cycle of these two sources through all fre- 
quencies taking about half an hour before being repeated. 
Observations were taken irregularly over the period 17-20 
June totalling approximately 36 hours of data. Full Stokes 
parameters were observed so that polarization maps could 
be made. 

Initial data reduction followed standard procedures for 
MERLIN data. After the flagging of corrupted data, the flux 
density scale was set relative to 3C 286. Subsequent cal- 
ibration was accomplished using the NRAO Astronomical 
Image Processing Software (aips). The phase-calibrator was 
found to be substantially resolved and was therefore itera- 
tively mapped and self-calibrated several times to produce 
an accurate model of the source structure. This was then 
used to derive antenna-based complex gain solutions which 
were then applied to the B0218+357 data. B0233+359 was 
also used to derive instrumental polarization corrections and 
3C 286 was used to calibrate the polarization position angle. 

The data from each frequency were next indivi dually 




Einstein Ring 



mapped an d self-calibrated using the difmap package ( 3hep 



herd 1997 ), solving for both phase and amplitude antenna 
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Figure 1. MERLIN single- frequency (4546 MHz) map made 
using the DIFMAP program. The data are uniformly weighted 
and contours are plotted at -0.2, 0.2, 0.4, 0.8, 1.6, 3.2, 6.4, 
12.8, 25.6 and 51.2 per cent of the peak brightness in the map 
(759 mjy beam - 1 ). The restoring beam of 51 X 47 mas is shown 
in the bottom-left corner. 



corrections. A map made from the 4546 MHz data is shown 
in Fig. [j]. Before the three frequencies could be combined 
into a single data set, it was necessary to correct the data 
for the steep spectral index of the Einstein ring (a = —0.6, 
where S v oc v a ). Simply scaling each frequency's data was 
not possible due to the compact cores having a significantly 
flatter spectral index than the ring (a — —0.2). This prob- 
lem was solved by subtracting the compact cores from the 
(u, v) data. This was achieved by (in aips) making a map of 
the calibrated data output from difmap and subtracting the 
Fourier transform of the CLEAN components correspond- 
ing to the positions of components A and B. The data were 
then remapped and the process repeated until all the core 
emission had been removed. The total flux density that re- 
mained in each subtracted map was then used to scale the 
(u, v) data sets to a consistent amplitude scale. This correc- 
tion was as expected, the higher frequencies having to be 
scaled upwards in amplitude relative to the lower. Finally, 
the CLEAN components removed from the central frequency 
were added back into the combined data set resulting in a 
'complete' B02184-357 data set. 

These data were then amplitude and phase self- 
calibrated (using the CLEAN components from an earlier 
4546 MHz map) and mapped. In order to remove spurious 
data, the brightest 1000 CLEAN components were vector- 
subtracted from the (u, v) data and points with a residual 
amplitude of >200 mjy were removed from the data. This 
cut-off was chosen using a "by-eye" estimate of the scatter 
in the data, but as these only constituted <0.7 per cent of 
the data the effect on the resulting map is minimal. The 
subtracted CLEAN components were then added back into 
the clipped data set, a map of which is shown in Fig. |^. In 
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Figure 2. MERLIN 5-GHz MFS map. Contours in this and sub- 
sequent maps are plotted at -3, 3, 6, 12, 24, 48, 96, 192, 384, 
768, 1536 and 3072 times the off-source rms noise in the map 
(102/iJy beam -1 ). The circular restoring beam of 50 X 50 mas is 
shown in the bottom-right corner. 



Figure 3. VLA 5-GHz map. The data are uniformly weighted 
and the bottom contour is equal to three times the off-source rms 
noise (60/iJy beam - 1 ). The restoring beam of 341 X 334 mas is 
shown in the bottom-right corner. 



this case the dirty map was made using MX and CLEANed 
using the task SDCLN. This task utilises the Steer-De wdney 



Ito (S DI) method of choosing CLEAN components ( Steer 
Dewdni ;y & Ito 1984 ) and is useful when a map contains ex- 



tended structure of similar brightness. This map represents 
a considerable improvement on the single- frequency maps 
(such as that shown in Fig. [j]), the dynamic range having 
increased from 1500:1 to 8000:1. The map includes some 
low-level brightness from the jet to the south of the Ein- 
stein ring. There are also obvious miscalibration artefacts 
(alternate positive and negative sidelobes) around compo- 
nent A. 



2.2 VLA data 

The VLA data were observed on the 8th August 1995 for a 
total of ~100 minutes, less than two months after the MER- 
LIN observations. Two separate frequency bands were used 
and these were centred at 4885 and 4835 MHz, each with 
a bandwidth of 50 MHz. 3C 48 was used as the primary 
flux calibrator and 3C 84 the point source phase-calibrator. 
Calibration followed standard procedures and proved to be 
straightforward with very little flagging of data necessary. 
B0218+357 was next iteratively mapped (using MX with uni- 
form weighting) and phase and amplitude self-calibrated to 
produce a map of the source which is shown in Fig. |§| This 
has an rms noise of 60 /iJy beam -1 , a figure that is substan- 
tially higher than the theoretical noise of ~16 ^tJy beam - . 
This though is unsurprising given that the dynamic range in 
the map is ~ 15 000:1. The map contains some likely spuri- 
ous features (particularly emission to the west of component 
A and the jet), but these are only present at the 3<r level, 
the lowest contour in the map. As the errors present in the 



VLA map are at such a low level of significance, they are 
unlikely to significantly degrade the quality of maps made 
from the combined MERLIN/VLA dataset. 



2.3 MERLIN/VLA combined data 

As the VLA data were observed in J2000 coordinates and 
MERLIN in B1950, the VLA data were transformed to 
B1950 coordinates. A map made from the transformed (u, v) 
data was essentially identical to that shown in Fig. |^. As a 
check on the accuracy of the coordinate transformation, the 
task maxfit was used to find the peak brightness and its 
position in both the MERLIN and the transformed VLA 
data sets, the results of which are shown in Table |l[ 

The VLA data were observed as close as possible in 
time to the MERLIN MFS data so as to keep to a minimum 
any complications in the combination of the two arrays aris- 
ing from source variability. This has been successful as the 
peak brightness in each map is very similar. The positions 
of the peaks are a little less consistent. Whilst in declination 
the peaks are coincident to within 1 mas, the right ascen- 
sion coordinates differ by 60 mas. Therefore the two maps 
are offset by about one MERLIN beam. This is not com- 
pletely surprising as no attempt was made to perform exact 
astrometry with the VLA observations. The poor resolution 
of the VLA map may also be a contributing factor. However, 
as the offset in position is not unduly large (and removable 
with self-calibration) and the flux scales broadly consistent 
the two data sets were simply combined with no scaling or 
removal of model components (A and B). Prior to this the 
weights of each visibility of each array were made approxi- 
mately the same so that each data set contributed roughly 
equally to the resultant image. 

All maps of the combined data set were made using 
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Table 1. Positions and brightnesses of the brightest points (component A) in the MERLIN (three-frequency) and VLA maps. 



Array Peak brightness (mjybeam 1 ) RA (hms) 



Dec. (° "') 



MERLIN 
VLA 



933.2 
930.3 



02 18 04.12603 +35 42 31.8348 
02 18 04.12998 +35 42 31.8335 



the MPS imaging task IMAGR. A value for the ROBUST pa- 
rameter (which allows a compromise to be made between 
the traditional natural and uniform weighting schemes) of 
— 1 was used in all maps which resulted in a beamsize of 
57 x 55 mas. The initial map of the combined data was very 
poor, but with several iterations of phase self-calibration 
subsequent maps were of much higher quality. In order to 
make the best possible map the data were also amplitude 
self-calibrated and corrected for baseline errors. This lat- 
ter step was particularly successful in removing the sidelobe 
structure around component A seen in Fig. ^. The final im- 
age is shown in Fig. ^ and has an rms noise of 82 /ijy beam^ 1 
and a dynamic range of 10 000:1. The dynamic range of a 
typical bright area of the ring is about 100:1. 



3 DISCUSSION 
3.1 Image fidelity 

The final image shown in Fig. ^ represents a marked im- 
provement on previous maps made of the Einstein ring in 
this lens system, combining the sensitivity of the VLA with 
the resolution of MERLIN. As the dynamic range of the map 
is greater and the aperture coverage so much better than in 
previous high-resolution maps, we can also expect there to 
have been a substantial improvement in the image fidelity 
(fractional on-source errors). Short of performing compli- 
cated and time-consuming simulations of the entire map- 
ping process on a model source it is difficult to calculate the 
image fidelity. Instead, in the following paragraphs we will 
consider several ways in which the theoretical image fidelity 
could be degraded and show that the magnitudes of these 
are negligible. 

The fundamental problem with MFS that must be over- 
come is that source brightness varies with frequency. We 
have compensated for this with the MERLIN data, as de- 
scribed in Section 2.1, by removing the flat-spectrum cores 
and scaling the remaining emission. In doing so we have 



assumed that the Einstein ring emission is described by a 
single value of a. As gravitational lensing is an achromatic 
process this is in general a reasonable approach, providing 
that the lensed source has a uniform spectral index. If this 
is not the case then errors will result from the spectral- 
index residuals. We believe that the assumption of uniform 
a holds fairly well in B0218+357 due to the fact that the 
area of jet imaged into the ring is small (of order 10 mas) 
and bec ause spectral-index g radients along radio jets are 
shallow (Bridle & Perley 1984). Furthermore, Conway et al. 



ignored when the dynamic range in the map is <1000:1. 
Our observations easily fulfill this criterion as the dynamic 
range of the ring emission is ~ 100:1 and the bandspread 
±7 per cent. The much brighter cores do not contribute to 
the spectral sidelobes as after their subtraction at all three 
frequencies, only those from the central frequency were sub- 
sequently returned to the combined data. As the average 
VLA and central MERLIN frequencies differ by less than 
1 per cent and the peak brightnesses in the three-frequency 
MERLIN map and VLA map were so similar, it is unlikely 
that major spectral errors could result from the addition of 
the VLA data. 

Another effect that will reduce the image fidelity is 
source flux density variability. This needs particular con- 
sideration with regards to B0218+357 as the radio core of 
the background source is variable, as it had to be for the 
time delay to be measured. Fortunately, the relatively low 
frequency of these observations means that any source vari- 
ability should be reduced compared with the rapid varia- 
tions seen at higher frequencies. VLA monitoring data at 
5{j] 8.4 and 15 GHz ( |Biggs et al. 1999| ) show that although 



(1990) have shown that for a knotty jet (of which the ring 
in B0218+357 could be considered an example with extreme 
curvature) observed with the MERLIN array with a band- 
spread of <25 per cent, the spectral errors can usually be 



highly variable at the highest frequency, the variations be- 
come much reduced in magnitude (by a factor of about one 
half) and have a longer timescale at 5 GHz. This theoretical 
argument is supported by the very similar peak flux densities 
in the MERLIN and VLA maps of this work despite their 
being observed approximately a month and a half apart. 



3.2 The new image 

Perhaps the most striking feature of Fig. ^ is the absence 
of emission in the middle of the ring lying very close to, 
and to the west of, component B. Related to this is the 
reduction in brightness that forms a valley running north- 
south through the ring and passing through this hole. A 
similarly-aligned feature was seen at low significance in the 
independently-mapped 5-GHz MERLIN image of Patnaik 
et al. (1993), lending confidence in the reliability of this fea- 
ture. The presence of this valley results in the ring taking 
on more the appearance of two arcs. Further support for the 
reduction in brightness close to component B comes from 
the 15-GHz VLA map shown in Biggs et al. (1999). The ra- 
dio jet can be seen emerging from the southern edge of the 
Einstein ring as the two protrusions at the bottom. These 
correspond well to similar features in the VLA 15-GHz radio 
map of Biggs et al. (1999). 

Another interesting area of the map is the bright spot 



t Although only the 8.4 and 15 GHz data featured in Biggs et al. 
(1999), a small amount of time per epoch was spent observing at 
5 GHz. 
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Figure 4. MERLIN/VLA 5-GHz map. The data are uniformly weighted with a ROBUST parameter of —1. The bottom contour is equal to 
three times the off-source rms noise (82/iJy beam - 1 ). Greyscales represent brightnesses between —1 and lOmJybeam - 1 . The restoring 
beam of 58 X 56 mas is shown in the bottom-right corner. 



(the third brightest component in the map) that lies out- 
side the Einstein ring, ~150 mas to the north of compo- 
nent A. This has not been identified before and is notice- 
able by its absence in the MERLIN 5 GHz map of Patnaik 
et al. (1993). Recently though, a feature at the same po- 
sition and with similar brightness has been seen in a com- 
bined MERLIN /European VLBI Network (EVN) image at 
1.7 GHz (A. R. Patnaik, private communication). In addi- 
tion to this, four-telescope VLBI data (Vermeulen et al., in 



preparation) shows HI absorption associated with the lens- 
ing galaxy against three components — A, B and a third 
component located close and to the north-east of A. For the 
purposes of this discussion we shall refer to this as compo- 
nent C. What this third component might be is currently 
undetermined, but given that its position is so close to com- 
ponent A it is unlikely that it is not a lensed image. Its 
counterpart (component D say) is most likely to be located 
close to component B, but if we assume that the relative 
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magnification between C and D is similar to that between 
A and B, component D will be difficult to detect. This is 
because it will be weakened by about a third compared to 
C and closer to component B than C is to A by a similar 
factor. 

It is the lensed emission in the Einstein ring that in- 
terests us most as it is this which will enable us to con- 
strain the lens model. One way of doing this is to identify 
lensed images of the same surface brightness as these are 
likely to be images of the same background source (as grav- 
itational lensing preserves surface b rightness). This is th e 



basis of the Ring Cycle algorithm ( Kochanek et al. 198!: ) 
that attempts to optimise the mass model so that, when 
back-projected into the source plane, the surface brightness 
dispersion within individual pixels is minimised. A natural 
consequence of this process is a reconstruction of the un- 
lensed source. However, it is clear from Fig. ^ that the two 
arcs of the Einstein ring contain substructure with differ- 
ent surface brightnesses. This is almost certainly due to the 
fact that the finite beam of the combined observations is 
not sufficient to resolve fully the brightness variations in the 
Einstein ring. I n this case the appropriate action is to use 
the LensClean ( Kochanek & Narayan 1992 ) algorithm. This 



represents a considerable improvement on the Ring cycle and 
uses the CLEAN algorithm familiar to radio astronomers to 
allow for resolution effects like those present in the combined 
MERLIN/VLA map of B0218+357. As with the Ring cycle, 
LensClean simultaneously finds the best-fit lens model and 
produces a map of the unlensed source structure. This work 
is currently underway and will feature in a future paper 
(Wucknitz et al., in preparation). 



4 CONCLUSIONS 

In this paper we have presented a new image of the gravi- 
tational lens system B0218+357 at a mean frequency close 
to 5 GHz. This has combined data from both the VLA and 
MERLIN, the latter in multi-frequency mode observing at 
three separate frequencies. The resulting data set has excel- 
lent (u, v) coverage and sensitivity, as well as resolution high 
enough (~ 50 mas) to easily resolve the Einstein ring. The 
new map is the best made of this system. 

The motivation for making this map was to improve 
the lens model for this system and hence to improve the es- 
t imate of Hp. This currently stands at 69^g kms" 1 Mpc -1 



(Biggs et al. 199£), but there are only limited constraints 



on the mass model available from the VLBI substructure 
of the compact images and the centre of the lensing galaxy 
is poorly determined. The next stage will be to exploit the 
new image for lensing constraints. The results of applying 
the Lens Clean algorithm will allow us to improve substan- 
tially the estimate of _ffo from the time delay. 



Network Contract ERBFMRXCT96-0034 'CERES'. MER- 
LIN is run by the University of Manchester as a National 
Facility on behalf of PPARC. The VLA is operated by the 
National Radio Astronomy Observatory which is a facility 
of the National Science Foundation operated under cooper- 
ative agreement by Associated Universities, Inc. 



REFERENCES 

Biggs A. D., Browne I. W. A., Hclbig P., Koopmans L. V. E., 
Wilkinson P. N., Perley R. A., 1999, MNRAS, 304, 349 

Bridle A. H., Perley R. A., 1984, ARA&A, 22, 319 

Browne I. W. A., Patnaik A. R., Walsh D., Wilkinson P. N., 1993, 
MNRAS, 263, L32 

Conway J. E., Cornwell T. J., Wilkinson P. N., 1990, MNRAS, 
246, 490 

Kochanek C. S., 1990, in Mellier Y., Fort B., Soucail G., 

eds, Gravitational Lensing, Lecture Notes in Physics 360. 

Springer- Verlag, Berlin, p. 244 
Kochanek C. S., Narayan R., 1992, ApJ, 401, 461 
Kochanek C. S., Blandford R. D., Lawrence C. R., Narayan R., 

1989, MNRAS, 238, 43 
Kormann R., Schneider P., Bartelmann M., 1994, A&A, 284, 285 
Lawrence C. R,., 1996, in Kochanek C. S., Hewitt J. N., eds, Proc. 

IAU Symp. 173, Astrophysical Applications of Gravitational 

Lenses. Kluwer, Dordrecht, p. 299 
Lehar et al., 2000, ApJ, accepted 

Patnaik A. R., Browne I. W. A., Wilkinson P. N., Wrobel J. M., 

1992, MNRAS, 254, 655 
Patnaik A. R., Browne I. W. A., King L. J., Muxlow T. W. B., 

Walsh D., Wilkinson P. N., 1993, MNRAS, 261, 435 
Patnaik A. R., Porcas R. W., Browne I. W. A., 1995, MNRAS, 

274, L5 

Refsdal S., 1964, MNRAS, 128, 307 

Sault R. J., Conway J. E., in Taylor G. B., Carilli C. L., Perley 
R. A., eds, ASP Conf. Ser. Vol. 180, Synthesis Imaging in 
Radio Astronomy II. Astron. Soc. Pac, San Francisco, p. 419 
Shepherd M. C, 1997, Astron. Data. Anal. Software Syst., 6, 77 
Steer D. G., Dewdney P. E., Ito M. R., 1984, A&A, 137, 159 



ACKNOWLEDGEMENTS 

ADB acknowledges the receipt of a PPARC studentship. 
We thank Olaf Wucknitz for fruitful discussions and the 
anonymous referee for several suggestions that significantly 
improved this paper. This research was supported in part 
by the European Commission TMR Programme, Research 



© 0000 RAS, MNRAS 000, 000-000 



